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We investigate the spherically symmetric gravitational collapse of an incoherent 
dust cloud by considering a LTB type space-time in third order Lovelock Gravity 
without cosmological constant, and give three families of LTB-like solutions which 
separately corresponding to hyperbolic, parabolic and elliptic. Notice that the con- 
tribution of high order curvature corrections have a profound influence on the nature 



PJ^ i of the singularity, and the global structure of space-time changes drastically from the 

analogous general relativistic case. Interestingly, the presence of high order Lovelock 

> 

f^ ' terms leads to the formation of massive, naked and time-like singularities in the 7D 

m ' 

t^^ ■ space-time, which is disallowed in general relativity. Moveover, we point out that 

(N 

f^^ . the naked singularities in the 7D case may be gravitational weak therefore may not 



be a serious threat to the cosmic censorship hypothesis, while the naked singulari- 
ties in the D > 8 inhomogeneous collapse violate the cosmic censorship hypothesis 



^ ■ seriously. 



I. INTRODUCTION 



One of the most intriguing predictions of Einstein's general relativity is the existence of 
black holes which are formed from gravitational collapse in the last stage of heavy stars' 
life or in high-density regions of the density perturbations in the early universe [l|. The 
gravitational collapse of an incoherent spherical dust cloud is described by the Einstein 
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equation Gab = £it,r)uaUh with the energy density e(t,r) and the time-hke velocity vector 
Ua- In a reference frame comoving with the collapsing matter, it can be proved that the 
metric of spherically symmetric space-time reads |2| 

ds^ = -de + A{t, rfdr^ + R{t, r)\de^ + sin^ ed(f)'^), 

(1) 

where r is the comoving radial coordinate and t the proper time of freely falling shells. As 
well known, the solution of this metric is the Lemaitre-Tolman-Bondi (LTB) solution J2|, 
which has been extensively studied not only in spherical collapse, but also in cosmology as 
well |3||4|. 

It is well known that the end state of CTavitational collapse is a singularity of space-time 
curvature with infinite density of matter[5| 6|. But, it is not known clearly whether such 
a singularity will be naked or covered by a event horizon. At the singularity, the laws of 
science and the ability to predict the future would break down. It led Roger Penrose to 

n n 

propose the (weak and strong) cosmic censorship hypothesis (CCH) |6|[7[, which to date 
remain unproven. The weak CCH asserts that there can be no singularity visible from 
future null infinity and light from singularity is completely blocked by the event horizon. 
It protects external observers from the consequences of the breakdown of predictability 
that occurs at the singularity, but dose nothing at all to observers who are cofalling with 
the collapsing massive body, while the strong CCH prohibits singularity's visibility by any 
observer. Despite almost 40 years of effort we still don't have a general proof of CCH |8|. 
On the contrary, it was shown that the seminal work of Oppenheimer and Snyder [l| is 
not a typical model and the central singularities formed in generic collapse are naked and 
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observable 

Up to now, many quantum theories to unify all fundamental interactions including gravity 
have been proposed, the most promising candidate among them is superstring/M-theory. In 
string theory, extra dimensions were promoted from an interesting curiosity to a theoretical 
necessity since superstring theory requires a ten-dimensional space-time to be consistent 



from the quantum point of view ll[[12|. It caused a renewed interest about the general 



relativity in more than 4 dimensions. Several LTB-li 

higher dimensions have been obtained in recent years 

years another renewed interest has grown in higher order gravity, which involves higher 



ike solutions to Einstein equations in 
higher dimensions have been obtained in recent years 13| [ij] • On the other hand, in recent 



derivative curvature terms. Among the higher curvature gravities, the most extensively 
studied theory is so-called Lovelock gravity 15|, which naturally emerged when we want 
to generalize Einstein's theory in higher dimension by keeping all characteristics of usual 
general relativity excepting the linear dependence of Riemann tensor. In addition, Lovelock 
terms naturally occur in the effective low-energy action of superstring theory 16| 17|. Hence, 
the Lovelock gravity provides a promising framework to study curvature corrections to the 
Einstein-Hilbert action. Since we are interested in what influence to the theory of gravity 
will be caused by Lovelock terms, for simplicity, we only consider pure Lovelock gravity 
and neglect the other terms in the low-energy expansion string theory. The Lagrangian of 
Lovelock theory is the sum of dimensionally extended Euler densities 15 | 
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where arbitrary constants «„ are the Lovelock coefficients, and £„ is the Euler density of a 
2 /c- dimensional manifold 
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Here the generalized delta function S^l...^^ ^^ totally antisymmetric in both sets of indices and 
R'^'^ab is t^^ Riemann tensor. Though the Lagrangian of Lovelock gravity consists of some 
higher derivative curvature terms, its field equations of motion contain the most symmetric 
conserved tensor with no more than two derivative of the metric. The Co is assumed to be 
identity and cq the cosmological constant. £i gives the usual Einstein-Hilbert action term. 
£2 is called Gauss-Bonnet term, it is a correction term of the action 18|. So far, the exact 



static and spherically symmetric black hole solutions in third order Lovelock gravity were 
first found in 19|, and the thermodynamics have been investigated in [20|][2l| 22|. 



It should be interesting to discuss the gravitational collapse in higher dimensional Love- 
lock gravity. The natural questions would be, for instance, how does the Lovelock terms 
affect final fate of collapse? What horizon structure will be formed? Whether solutions leads 
to naked singularities? Whether CCH can be hold? Recently, Maeda considered the spheri- 
cally symmetric gravitational collapse of a inhomogeneous dust with the D > 5-dimensional 
action including the Gauss- Bonnet term. He discussed its effects on the final fate of grav- 
itational collapse without finding the explicit form of the solution 23|]. Then, Jhingan and 
Ghosh considered the 5D action with the Gauss-Bonnet terms for gravity and give a exact 



model of the gravitational collapse of a inhomogeneous dust [2J][25|]. It's interesting to in- 
vestigate the effect of higher order terms of Lovelock gravity in gravitational collapse, we 
explore the gravitational collapse in comoving coordinates, and seek LTB-like solutions in 
the third order Lovelock gravity. The dimensions of space-time concerned is D > 7 because 
the third order term yields nontrivial effect in dimensions greater than or equal to 7 26| 27 1. 
Using our solutions, we discuss the formation of singularities, and analyze the nature of 
them. In particular, we consider that whether such singularities would be hidden or be 
visible to outside observers. Since the CCH would be violate by naked gravitational strong 
singularities, we investigate the gravitational strength of naked singularities, and discuss 
that whether the naked singularities in third order theory is a serious threat to CCH. 

This paper is organized as follows. In Sec [TTl for the D > 7-dimensional space-time, we 
give the field equations in third order Lovelock gravity without a cosmo logical constant, 
and derive the LTB-like solutions. In Sec llllt we investigate the final fate of the spherically 
symmetric gravitational collapse of a dust cloud. The subject of Sec|IV]is to analyze horizons 
(both apparent horizons and event horizons) and trapped surfaces, and explore that whether 
the singularity formed by gravitational collapse is hidden or visible to outside observers. The 
strength of the singularity is demonstrated in Sec |Vl Sec I VI I is devoted to conclusions and 
discussions. 

Throughout this paper we use units such that 8nG = c^ = 1. 

II. LTB-LIKE SOLUTIONS IN THIRD ORDER LOVELOCK GRAVITY 

The action of third order Lovelock gravity in the presence of matter field can be written 



as 



S = d^Xy/^{R + 02-^2 + asA) + Smatter, (3) 

where 02 and as are coupling constants of the second order (Gauss- Bonnet) and the third 
order terms, respectively. In the low-energy limit of the heterotic string theory, a is regard 
as the inverse string tension and positive define [28]. Hence we restrict ourselves to the 
case a > in this paper. For future simplicity, we take coefficients 02 = id-3)'(d-4) ^^^ 



as = 72C^ — • The Gauss-Bonnet term £2 is 
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and the third order terms of Lovelock Lagrangian is of the form 
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Varying the action Eq. Q, we obtain the equation of gravitation field 

Gi:> + a^Gf, + a,Gf, = T,^, (4) 

where 

G'(^) = R - -Ra 
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The solution we search is collapse of a spherically symmetric dust in D > 7 space-time. 
Following LTB solutions, we assume that the system consists of a freely falling perfect fluid, 
it requires that the mean free path between collisions is small compared with the scale of 
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lengths used by observer. The energy-momentuni tensor for the perfect fluid in comoving 
coordinates is 

T^^ = e{t,r)u^u^, (5) 

where it^ = 5* is the velocity vector field. We assume that the energy density e{t, r) on the 
initial surface is smooth, that is, it is can be extended to a C°° function on the entire real 
line. The metric in comoving coordinates is written in the form 

ds'^ = -dt^ + A{t, rfdr^ + R{t, rfdQl_2, (6) 

where r is the comoving radial cooedinate, and t is the proper time of freely falling shells. 
Plugging the metric Eq. IQ into Eq. (jl]), we have 

Gl = '^^{AE! - AR!)[[5{A^F? + A^ - H^f 

+ 2aA^R^{A^R^ + A^ - R'^) + A^R^] = 0, (7) 

where an over-dot and prime denote the partial derivative with respect to t and r, respec- 
tively. Based on Eq. (jTj), we arrive at two families of solutions which satisfy 

= ii?' - AR!, (8) 

and 

= ^{A^ie + A^- R^ 

+ 2aA^R\A^R? + A^ - R'^) + A'R\ (9) 

respectively. The second equation involves the Lovelock coupling constants, and leads to a 
trivial solution if a — )■ 0, /3 — > 0. We will neglect it since we want to explore the Lovelock 
corrections to Einstein gravity. The solution of Eq. ([8]) reads 

with an arbitrary function W{r). 

As it is possible to make an arbitrary re-labeling of spherical dust shell by r — )■ g{r), we 
fix the labeling by requiring that, on the hypersurface t = 0, r coincides with the area radius 

i?(0,r)=r. (11) 



Apparently, every t = const and r = const slice of the space-time is a sphere of radius R. 
The radius R can be given an absolute significance by following interpretation. Considering 
two particles a and b distribution along the radial direction in comoving coordinates, the 
space distance between them at the coordinate time t can be obtained on the hypersurface 
t = const as 



PO PD 

lab = / dl = ^Jhijdx^dx^ 

J a J a 



dr= / , (12 

here hab is the induced metric on the hypersurface, K{r) = iy(r)^ — 1 and we assume 
K{r) > — 1. In last step the fact dR = R'dr + Rdt = R'dr at a constant t is used. Thus, 
the line element of actual space distance along radial direction is ,^^ . The signature 
of K{r) corresponds to three types of solutions, namely hyperbolic, parabolic and elliptic, 
respectively. K{r) = is the marginally bound case in which the metric takes the form of 
Minkowski metric on the hypersurface t = 0. 

In comoving coordinates, the equations of momentum conservation {T^)-^^ = are auto- 
matically satisfied, and the t-component reads 

n = _d^_.{R;^.{D-2){R^ 
dt ^ 2i?'2 2/?2 ^ 

= -|M^-^i?'), (13) 

which gives the solution 

Thus, the mass function is defined as 



Mir) = .^— j e{t,T)R''-'dR 
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e(0,r)r^-2rfr, (15) 



D-2 

it is positive and increases with increasing r for the nonnegativity of energy density e. 
Based on Eq. ( ITO|) and K{r) = W{ry — 1, the tt component of equations is given by 

2-D 



^* " QR^-^R' 



R^-^((3{R^ - Kf + 3a(i?2 - KfR^ 



+ 3{if-K)R^) =-e{t,r). (16) 



Substituting Eq. ( !T^ and Eq. (fT5|) into Eq. (fT6|) . we get the real solution that 
where 



^^ /22 e + n.i ,e-n,i , 



(17) 



n = 



/32(/32p2^ - 4:a^pD + Qal3pD - Sa^ + 4/3) 

with p£) = MR^~^ . This equation governs the time evolution of R in D-diniensional third 
Lovelock gravity. For gravitational collapse, the solution of -R(t, r) takes 



r=-JK+ 



T 



+ I — z — ) -a 



(18) 



It is straightforward to check that other field equations are automatically satisfied when 
Eq. (fTTl) is satisfied. 

For arbitrary initial data of energy density e(0,r), the Eq. (IT8l) completely specify the 
dynamical evolution of collapsing dust shells. In the general relativistic limit a — > 0,/3 — )■ 
with marginally bound case, Eq. (fT8|) can be integrated to yield 

2V3R- 



h{r) - t{r) 



2 



(D-l)v/M(0' 
where tb{r) is an function of integration, and can be formulated as 

2y3rV 



(19) 



tb{r) 

Consequently, the function R(t, r) is 

R{t, r) = r 



(D-l)v/MW' 



1 + ^^ i_, t^/m^) 



2 
D-1 



(20) 



2y3r^ 

It is similar to the form of R{t, r) in LTB and LTB-like-solutions |2||l4|. 

We can consider that our LTB-like solution is attached at the boundary of the dust cloud, 
which is represented by a finite constant comoving radius r = tq > 0, to the outside vacuum 
region. The outside vacuum region is represented by the solution whose metric is 



^2 



ds' = -F{r)dT' + ^ + r'dQl_„ 
V ^ F(r) ^ ^ 



where F(r) takes form of 



where 
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e = -2a^ + 3af3 + f3^pD, 




n = ^/32(/32p2^ - Aa^pD + QaPpD - 3a^ + 4/3) 
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v'l + '<'('■'.*, 



with p£) 

m 
R 

we can prove that the LTB-like solution attached at a finite constant comoving radius r = tq, 
where we represent this hypersurface as S, to the outside vacuum solution smoothly. 
Proof: As seen from inside of S, the metric on E is obtained by 

dsl = -df + i?|dfi^_2. 

As seen from outside of E, it is 

dsl = -iFiR^)Ti - -^)dt' + Rld^l.,. 

It can be checked that the induced metric is the same on both sides of the hypersurface E 
with the definition Eq. fl22|) . It implies 

F{R^)Tl - -^ = 1. (23) 

We can define a function C, by 



C=dF{R^) + Rl = F{Rj,)Tj,. 



As seen from inside, the nonzero components of the extrinsic curvature K,^ of E are calculated 
as /C* = and /Cj = ^^^^^. As seen from outside of E, we find /C* = X(F(i?s)T| - 
p,^ s )~^ and K,\ = -^. With the definition Eq. ( |22l) . /C^ is the same on both sides of E. It 
implies 

C = y/l + K{ro). (24) 
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Finally, we combine Eq. ( 123|) and Eq. (^^ giving the equation of motion for the hypersurface 

S as 

Rl = l-F{Rj:) + K{ro). 

It takes the same form as Eq. fl29|) with the definition Eq. fl22l) . Thus, the LTB-like solution 
attached at a finite constant comoving radius r = ro to the outside vacuum solution smoothly. 

D 

III. SHELL FOCUSING SINGULARITY 



In this section, we consider the final fate of the gravitational collapse in D > 7 space-time. 
As pointed earlier, in general relativity, a collapse leads to a singularity, and the conjecture 
that such a singularity must be covered by an event horizon is the weak CCH. There are 
two kinds of singularities: shell crossing singularity and shell focusing singularity which is 
defined by i?' = and R = 0, respectively. The characteristic of a singularity in space-time 
manifold is the divergence of the Riemann tensor and the energy density 6| . In our case, 
the Kretschmann invariant scalar /C = R^^^rR'^"'^^ for the metric Eq. i^ is 

2R'^ A{D-2)R^ A{D-2)R^R'^ 

It can be certified that the Kretschmann scalar is finite on the initial data surface. According 
to Eq. (TT5l) , the energy density of fluid dust sphere is 

sit r) - ^^ - ^)^^ (26) 

Clearly the Kretschmann scalar and the energy density diverge when i?' = and i? = 0. 
Hence, we have both shell crossing and shell focusing singularities. 

Shell crossing singularities can be naked, but they are inessential. Although the mass 
density and curvature invariants blow up there, the metric is in fact continuous. This is seen 
from that the Riemann curvature tensor and the energy momentum tensor are well defined, 
if we replace the coordinates (t, r) by {t,R). On the other hand, shell focusing singularities 
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Is 



are considered to be the genuine singularities in space-time manifold [9|][29|]. Henceforth, we 
only concern the shell focusing singularity here. 

In third order Lovelock gravity, equations and solutions of gravitational collapse are quite 
different from counterparts in general relativity. Hence it is necessary to investigate whether 
the evolution of collapsing dust cloud leads to the formation of the shell crossing singularity. 
We assume each shell satisfies R{0,r) < initially. If i? = is satisfied for some r after 
the initial moment, this shell ceases to collapse and then bounces {R > 0), and the shell 
focusing singularity would not be formed. Indeed, it can be proven that in the case /3 > 0, 
R^ — K > aX the initial moment, the shell focusing singularity will be inevitably formed. 

Proof: From Eq. ( IT6l) . we find that 

d{R^) _ (7 - D)MR^-^ - Qae'^R - 12iR^ 
dR 3/3£2 + QaiR^ + 3i?4 ' ^ ^ 

where (. = R? — K. Clearly, -W^ < if /3 > and £ > 0. Here we have used the condition 
a > which is mentioned in the last section. If d, is non-negative at the initial surface, it will 
be always non-negative as increasing R^. Hence, \R\ increase as decreasing R therefore the 
collapsing of the shell r = const from the finite initial data -R(0, r) = r is accelerated. Thus, 
R will vanish in finite proper time for the comoving observer, the shell focusing singularity 
is inevitably formed. D 

Physically, the condition Z > is satisfied for /3 > a^. For < /3 < a^, one should take in 
mind Eq.( ll7l) with a real solution, which gives a lower limit of po- Within this region, the 
numerical figure show the condition valid, but it is hard to analytically prove it. Moreover, 
one can show the condition / > valid forever if it is satisfied at the beginning. Thus, the 
final fate of a freely falling fluid sphere, which have initial density £:(0, r) and zero pressure, 
is a state of infinite energy density and curvature. In the rest of this paper we will only 
consider the case that the shell focusing singularity could be formed, and assume /3 > 0. 

As demonstrate above, our solution can be attached to the outside vacuum solution at 
the boundary of the dust cloud. The outgoing property of the central singularity depends 
on the dominant term of F{R) in vacuum solution for /? — )■ [34!]. If the metric function 
F < 0, then the tortoise coordinate defined by 

R*= I F~^dR 



has finite negative value, the singularity is space-like. If F > 0, the singularity is time-like. 
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If F = 0, the singularity is null since |-R*| — )■ oo. From Eq. (12T!) . we find that the metric 
function behaves around the central singularity as 

where r^ is the boundary of the collapsing dust cloud. Since we assume /3 > 0, if -D > 7, we 
have F{R) < hence the singularity is space-like. \i D = 1 and M{rij) > /3, the singularity 
is also space-like, li D = 7 and M{rf,) = 13, it is null, li D = 7 and M{ri,) < 13, it is time-like. 
This conclusion is shown in the Penrose diagram FIG.l. The event horizon of such vacuum 
solution would be discussed in the next section. 

If the energy density is independent on r, namely the energy density and the space 
are homogeneous, the time of the formation of the singularity tsF is a constant, as in the 
general relativity case. A homogeneous space that is isotropic about some point is maximally 
symmetric and the curvature in such a space is a constant. In homogeneous case, the metric 



Eq. ([6]) takes the form as the Robertson- Walker metric 30 1 



ds' = -de + a{tf{^^^ + r'dnl_,), (28) 

where /c is a constant. In order to satisfy Eq. f lTT|) . the radial coordinate r can be normalized 
so that a(0) = 1. When the shell hits the shell focusing singularity, the time is completely 
determined by a{t) = 0. Considering metric Eq. (128|) and M{r) = jjjzSijZii^i^Y^^^ i we 
get 



/, a(t)^,e + n,i ,e-n,i 

«W = ~\lk + ^[(^— )' + (^— )' - "]' (29) 



where 



e = -2a^ + 3a/3 + /J^e, 

n = y/P^{P^e'^ - Aa^e + Qa^e - Sa^ + 4/3) 

with e = 6e{0)a{ty^^ / {D — 1){D — 2). This equation does not contain the variable r, that 
is, a{t) is a function of t independent on r. It implies that every shell will collapse into the 
shell focusing singularity at the same time. Obviously, this conclusion is hold in general 



relativity fl|, Gauss-Bonnet gravity 2^ and third order Lovelock case. 

From the geometric perspective, it can be proven that in the homogeneous case the shell 
crossing singularity is ingoing null when D = 7 and is ingoing space-like when D > 7, as 
shown in FIG.l. 
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Proof: Considering Eq. (129|) we have 

6 



e(0) = [i3{d^-ky + 3a{d^-k)a 



{D-l)iD-2) 

+ 3a^](d2-A;)a^-'. (30) 

Since the energy density £:(0) have a finite non- negative value, we find that the factor a 

behaves as c{t — tsF)^~^ near the singularity t = tsF, where c is a coefficient. We take the 
line element of the FRW solutions with the factor a = c{t — tsp)^^ to the conformally flat 
form as 

ds^ = a\t{i, f))h\r{f)){-de + df^ + fHnl_^), (31) 



where 



dt 
dt 



a{t)h{ry b{r) 



hir) = exp(ln Irl + In I — 1-\/^ + ^|)- 

The range of t for t G {—oo, tsp) is (— oo, +oo) when D = 7 and is (— oo,to) for D > 7, 
where to is a constant. Thus, the shell focusing singularity is ingoing null for D = 7, while 
it is ingoing space-like for D > 7. n 

IV. HORIZON AND TRAPPED SURFACE 



An important construction in general relativity is that of the trapped surface, which is 
indispensable in proving null-geodesies incompleteness in context of gravitational collapse 
61]. In general relativity, the trapped surface is defined as a C^ closed space-like two-surface 
that two families of null geodesies orthogonal to this surface are converging. Physically, it 
captures the notion of trapping by implying that if two-surface T{t,r) {t,r = const) is a 
trapped surface, then its entire future development lies behind a horizon. The apparent 



horizon is the outermost marginally trapped surface for the outgoing null geodesies 32| . 
therefore the trapped surface could not be formed during coHapse without the occurrence 
of the apparent horizon. The main advantage of working with the apparent horizon is that 
it is local in time and can be located at a given space-like hypersurface. Instead, the event 
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(a)!) > 8 



{c)D = 7, M{n) = /3 



{d)D = 7, A/(rfc) < (3 



(b) 
i? = 7, Af (rfc) > (3 

FIG. 1: Penrose diagram of the homogeneous collapse of a spherically symmetric dust cloud in third 
order Lovelock gravity. Zigzag lines represent the shell focusing singularities, /+(~) corresponds to 
the future (past) null infinity, BEH and CH stand for the black hole event horizon and the Cauchy 
horizon, respectively. 

horizon coincide in case of static or stationary space-time, it is non-local. Moreover, in the 
vacuum region the apparent horizon coincides with the event horizon of the vacuum solution 
9| . That is, without the presence of apparent horizons there is no event horizon. 

Demanding the presence of the trapped surface in our spherically symmetric case implies 



g^'R^R, = -R' + — <0 



(32) 



and the condition for the existence of the apparent horizon with outward normals null is 



g^''R^,R,u = -R" + 






0. 



Using Eq. (fTOj) . the apparent horizon condition becomes 

R^ -K = l, 
Combing Eq. ( 1T6|) and Eq. ( l34l) . we obtain 

RAH(tAH{r),r) = 
where 



(33) 



(34) 



-3a + ^ 



6 



(35) 



^ = ^12M{r)RAH{tAH{r),ry-'' + 9a^ - I2f5. 



(36) 
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Clearly, coupling constants a and /3 produces a change in the location of apparent horizons. 
Such a change could have a signification effect in the dynamical evolution of these horizons. 
It has been shown that in the 5D Gauss-Bonnet gravity case, positive a forbids apparent 
horizon from reaching the coordinate center thereby making the singularity massive and 



eternally visible 24], which is forbidden in the corresponding general relativistic scenario 



Is 



and D > 6 Gauss-Bonnet gravity J9|] [23| [33| . In our case, positive a leads to noncentral naked 
singularities when D = 1. From Eq. fl36l) . we find the condition that the apparent horizon 
is earlier than singular shell {RAHitAHij)^^) > 0) for positive a in 7D case is M{r) > (3. 
Oppositely, Eq. ( l36l) could not be satisfied before the formation of the singularity in the 7D 
space-time in the case M{r) < /3. Furthermore, the apparent horizon touches the singularity 
when M{r) = /3. Thus, there is no shell can reach the apparent horizon if 

M{n) < /3, (37) 

with Tb the boundary of the dust cloud. It implies if the mass function M{rb) takes sufficiently 
small value, the shell focusing singularity is eternally visible from infinity during collapse, 
and leave open even the weak form of the CCH. 

The condition of the formation of the eternally visible shell focusing singularity in 7D 
case is completely determined by the initial data of the energy density. With the help of 
Eq. ( 1T5|) . we can find the initial data of the energy density that condition Eq. ( 1371) requires. 
For homogeneous case, such initial data satisfies 

^(M < p. (38) 

5 

We can consider a more realistic model that £:o[l ~ (~)"]) which is a density profile where 
energy density decreases as an observer move away from the center, as is expected inside a 
star. This form of initial data leads to 

5(n + 6) 

On the other hand, it is clear to see that in the D > 8 space-time, apparent horizons 
lies earlier than the singularity unless r = 0. Thus, the solution forbids the formation of 
the naked singularity except the singular point at the coordinate center, as same as in the 
general relativity and the Gauss-Bonnet gravity. As pointed by Christodoulou, the center 
singularity can be naked for suitable initial data of e(0,r) in the inhomogeneous case |9|. 
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The physical picture captured here is that the coordinate center hits the singularity so early 
that other shells have not reached apparent horizons at this moment, hence the light from 
the singular coordinate center can escape to the infinity. 

In Gauss-Bonnet gravity, the formation of the eternally visible shell focusing singularity 
is forbidden in D > 6 space-time, and is permitted when M < a and a > in the 5D case. 
Comparing results in second and third order Lovelock gravity, one can conjectures: for n 
order Lovelock gravity, such a naked singularity does not exist in D > 2n + 2 space-time, 
and would be formed when M < a„ in D = 2n + 1 case, where a„ is the coefficient of the 
highest order Lovelock term and is positive defined. 

Now we consider the event horizon of the vacuum solution which describes the final fate 
of the collapse. Based on Eq. ( 12T|) . we find the metric function F{R) satisfies 



2d2 



M{n) = R""'' I3{1 ~ F{R)y + 3a{l - F{R)yR 
+ 3{1 -F{R))R^ 

The event horizon occurs when F{Reh) = 0, thus we obtain 

M{n) = i?f^^(/3 + 3aRlj, + 3R%j,) 

This equation leads to that 



R 



EH 



\ 



3a + \ 12M{n)R'jfjy + 9a2 _ 12/3 



6 

Comparing this equation to Eq. ( 1361) . it is obvious to see the condition for the existence of the 
event horizon is the same to the condition for the presence of the apparent horizon, which is: 
it exist in the D > 8 case, and in the 7D case with M{rb) > (3. If M{rh) = /3 in 7D space-time, 
the event horizon touches the singularity at the center of the vacuum solution coordinate. 
This result is also shown in FIG.l. Thus, the D = 7 and M{rb) < /3 case indicates the 
existence of naked singularities which are eternally visible from infinity, and violates even 



the weak form of CCH. Such naked singularities are forbidden in standard general re 



ativity. 



and is allowed in collapse in 5D second order Lovelock gravity with positive a [23|. Such 
conclusion coincides with the result which is obtained by investigating the apparent horizon 
and the outgoing property of the singularity. 
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V. STRENGTH OF SINGULARITY 



In the case of formation of a naked singularity in gravitational collapse, one of the most 
significant aspects is the strength of such a singularity in terms of the behavior of the 



gravitational tidal forces in its vicinity [35||. The importance of the singularity strength 
lies in the fact that even if a naked singularity occurs, if it is gravitationally weak in some 
suitable sense, it may not have any physical implications and it may perhaps be removable 



by extending the space-time through the same [36|]. The gravitational strength of a space- 
time singularity is characterized in terms of the behavior of the linearly independent Jacobi 
fields along the time-like or null geodesies which terminate at the singularity. In particular, a 
causal geodesic 7(s), incomplete at the afiine parameter value s = Sq, is said to terminate in a 
strong curvature singularity at s = Sq, if the volume three-form V{s) = Zi{s) /\ ^'2(5) /\ Z-^^s) 
defined as a two-form in the case of a null geodesic ) vanishes in the limit as (s — ?■ sq) for all 
linearly independent vorticity free Jacobi fields Zi(s), ^2(5), ^3(5) along 7(5). A sufficient 
condition for a strong curvature singularity in 4D space-time is that in the limit of approach 

nn 

to the singularity, we must have along at least one causal geodesic 7(5) [29[[37[, 

lim (s - softfj = lim (s - sofRabVV'' > 0, (40) 

where V^ is the tangent vector to the geodesic. Essentially, the idea captured here is that in 
the limit of approach to such a singularity, the physical objects get crushed to a zero size, and 
so the idea of extension of space-time through it would not make sense, characterizing this to 
be a genuine space-time singularity. This condition has been applied to higher dimensional 



Here we follow them to assume that such 



case by Ghosh, Beesham and Jhingan 
condition can be applied to our case. 

We consider radial time-like causal geodesies U^ = ^ with the marginally bound case 
K = 0, here the affine parameter r is the proper time along particle trajectories. According 
to such definition, U^ satisfies If^U^ = —1, that is, 

- (^Y + R'y^f = -I, (41) 

dr dr 

Using the geodesic equation, we have 



iR")A-^y = 0. (42) 



dr^ 2' '''' dr 
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Substituting Eq. ( BTj) into Eq. (H2|) . we obtain that radial time- like geodesies must satisfy 

f4((^f-l)^0. (43) 

The simplest solution is the worldline of a freely falling particle, which is f/^ = ^ = 5f . In 
terms of proper time we can describe it as 

tsF{r) -t = To-T. (44) 

Eq. (13 7p shows that the naked shell focusing singularity in 7D space-time occurs when the 
region of the collapsing dust is sufficiently small. Moreover, as mentioned earlier, in Z^ > 8 
space-time only the central singularity can be naked. Hence, our purpose is to discuss the 
strength of the central singularity. We consider the expansion oi e{r) near r = 

+ 00 

£{r) = ^£nr'' ^eo, (45) 

ra=0 

it specify that R{t,r) behaves as Ro{t)r (homogeneous case) near the coordinate center, 
where -Ro(^) is a function of t and vanishes at t = tsF- Thus, we get 

^ = R^,U^U' = _i£_l)^. (46) 



Clearly, 



thus we have 



hm ^^ = hm ^°(^) - ^°(^^^) = lini R,(tl 

T^To T — Tq t"-s>to t — tgp T"-S>TO 



J^2 

lim (r - Tof = lim -^, 



^ jV=lim-(D-l)^° 



hm (r - ToY^ = hm -{D - 1)^. (47) 



From Eq. flT7|) . we obtain that 



limi?^= ii^( 0)3i?^— , 

T->-To T-S>ro p 

7 - D Mr, 1 luS 

limi?o= lim^-— (-^)3i?,3 (48) 



M{r) 

CU. CXD IV-LQ — mil 

Eq. (l47|l we have 



where Mq is defined as Mq = lim -^tr therefore has finite value. Substituting Eq. ( 148|) into 



(D-l)(D-7) 



lim (r - ro)'^ = ^ '-^ '-, (49) 
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the strong curvature condition Eq. (HOj) is not satisfied on the singularity near r = in the 
7D space-time, while it is satisfied in the D > 8 space-time. 

Hence, the central shell focusing singularity is gravitational strong in the D > 8 space- 
time and may be gravitational weak in the 7D space-time since we chosen a special solution 
of Eq. (I43l) . The difference between two cases is that the singularity can be eternally visible 
in the 7D case while it is forbidden in the D > 8 case. That is, the singularity is gravitational 
strong when the solution represents the black hole formation for arbitrary initial data. If a 
naked singularity is gravitational weak, it may not have any significant physical consequences 
so may not be a serious threat to the CCH. Thus, the naked singularity in the D > 8 case 
violates the CCH seriously, as in the general relativity 9l]J29l], while the naked singularity 
may not be regarded as an essential counter example to the CCH in the 7D case. 

Nevertheless, such a central singularity in the 7D case is gravitational strong in the general 
relativistic and second order Lovelock gravity. In the general relativistic limit a — )■ 0, /3 — )■ 0, 
Eq. ( I40l) takes the form 

hm (r - To) ^ = hm (r - Tq) [ ^^.d-i ]' 



with 



'T^ro T-S>To 6/?o 



3i?n^-^ 



2 _ 1- ^^0 



lim (r — To) = lim 

T-5-T0 T-S>ro Mi 





Thus, 

.im(r-.„ft>^ (^-^f-^) >0. (50) 

the strong curvature condition is satisfied in D > 3 space-time. In second order Lovelock 
gravity limit /3 — ;■ 0, we find that 

.■n.(r-^)V^ '^-^f-^) >0. (51) 

T-^TQ 4 

the strong curvature condition is satisfied in Z^ > 5 space-time. These results indicate 
that the Lovelock interaction weaken the strength of the singularity. Moreover, a generic 
expression for n order Lovelock gravity could be 

, , ,, (D-l)(D-2n-l) , , 

lim (r - T.fij = ^ '-^ i > 0. 52 

T-s>To An 

It shows that the strong curvature condition is satisfied in D > 2n + 1 space-time. 
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VI. CONCLUSIONS AND DISCUSSIONS 



In this paper, we have investigated the gravitational collapse of a spherically symmetric 
dust cloud in D > 7 space-time in third order Lovelock gravity without a cosmological 
constant. From field equations given by the third order Lovelock action, we discussed the 
solutions of three families K{r) > 0, K{r) < and K{r) = 0. We discussed the final fate 
of the collapse, and gave a condition for the formation of the shell focusing singularity. We 
also analysed the global structure of the space-time and gave the Penrose diagram for the 
homogeneous case. It turns out that, the contribution of high order curvature corrections 
has a profound infiuence on the nature of the singularity, and the whole physical picture of 
the gravitational collapse changes drastically. High order Lovelock terms alters the course 
of collapse and the time of formation of singularities, modifies apparent horizon formation 
and the location of apparent horizons, and changes the strength of singularities. 

The most attractive consequence is that an massive naked shell focusing singularity is 
inevitably formed in 7D space-time, which is quite different from that in the general rela- 
tivity and in the Gauss-Bonnet gravity. However, as we shown, the strength of the naked 
singularity in 7D case is weaker than that in the general relativistic limit, therefore this 
may not be a serious threat to the CCH. On the other hand, unlike the 7D case, there is 
a serious threat to the CCH caused by the naked, gravitational strong singularity in the 
D > 8 inhomogeneous collapse, thus the CCH which is violated in general relativity could 
not be protected in Lovelock gravity. 

When revising our paper in lines with reviewer's suggestions, we have received a paper 



by Ohashi, shiromizu and Jhingan 39| discussing the gravitational collapse in the Lovelock 
theory with arbitrary order which partially covers the result of present paper. In fact, 
they did not investigate the ingoing and outgoing properties of the singularity, the Penrose 
diagram, and whether the condition of strong singularity is satisfied. One may study such 
properties in a general Lovelock gravity similarly. 
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